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Effects of electron correlation on the electronic structure and magnetic properties of the Gd(OOOl) 
surface are investigated using of the full-potential linearized augmented plane wave implementation 
of correlated band theory ("LDA+U"). The use of LDA+U instead of LDA (local density approx- 
imation) total energy calculations produces the correct ferromagnetic ground state for both bulk 
Gd and the Gd surface. Surface strain relaxation leads to an 90 % enhancement of the interlayer 
surface-to-bulk effective exchange coupling. Application of a Landau-Ginzburg type theory yields 
a 30 % enhancement of the Curie temperature at the surface, in very good agreement with the 
experiment. 



Many of the magnetic properties of Gd metal are well 
understood Q. The half-filled 4f-shell = f,L = 0) 
of Gd leads to a formation of a well localized spin-only 
magnetic moment. These localized spin moments couple 
through an RKKY-type exchange interaction to form a 
ferromagnetic (FM) Heisenberg system with a bulk Curie 
temperature [T^!) of 293 K FM order polarizes the 
conduction electrons and leads to a total magnetic mo- 
ment of 7.63 /is / Gd atom 

However, in spite of relatively simple bulk magnetic be- 
havior, the magnetism of the Gd surface is rather unusual 
1^. The results of different spectroscopic measurements 
suggest a significant enhancement of the surface Curie 
temperature (T*) for Gd(OOOl). Gd is thus one of only 
three ferromagnets (including Tb and FeNis) for which 
such an increased has been observed [Q . After a first 
observation of this effect in Gd by Rau et al. Q], Weller 
et al. verified it for 400 A thick Gd(OOOl) films grown on 
a W(llO) substrate, by comparing spin-polarized low en- 
ergy electron diffraction (LEED) and the magneto-optic 
Kerr effect measurements and further suggested a 
possible anti-ferromagnetic (AFM) alignment of the sur- 
face layer(s) with respect to the bulk FM Gd. Further in- 
vestigations with spin-polarized valence and core photoe- 
mission spectroscopy (PES) |^ did not confirm an exis- 
tence of this surface AFM coupling. Instead, the in-plane 
component of surface layer magnetization was observed 
to be parallel to the bulk |^|| , althrough the possibility 
of canted or mixed in-plane and out-of-plane surface mag- 
netic ordering was suggested. Very recent spin-polarized 
photoelectron diffraction experiments |9| for bulklike ~ 
300 A thick epitaxial Gd/W(110) films clearly indicate 
temperature dependent core level spin asymmetries well 
above the bulk T^, also suggesting surface enhancement 
of of as much as 85 K. However, to date, there has 
been no quantitative theoretical explanation for this en- 
hanced T^. 

The aim of this paper is to show that first principles 
calculations which account for both electronic correla- 
tions of the 4f electrons and the relaxation of the sur- 
face atomic positions can provide such a quantitative de- 
scription of the electronic and magnetic structure of the 
Gd(OOOl) surface. Using the results of total energy cal- 



culations and the Landau-Ginzburg model we show that 
there is an increase of at the Gd surface due to an en- 
hanced surface-to-bulk effective exchange coupling that 
is in turn caused by the surface structure relaxation. We 
emphasize the substantial role of electron correlations for 
the Gd f-electrons to obtain a correct FM ground state 
for both bulk Gd and the Gd surface. We also emphasize 
the role of structure relaxation as a "driving force" for 
the Tc enhancement. 

1. Difficulties of local (spin) density approximation 
(LDA) and generalized gradient approximation (GGA) 
and limitations of a core treatment for Gd 4f electrons. 
Since the pioneering work of Dimmock and Freeman |ic| ] 
there have been other attempts to describe the 4f states of 
Gd in terms of the localized "4f-core" electron model jl^] 
(in which the 4f states of Gd are treated as a part of fully 
localized atomic core). Singh performed a detailed 
analysis of the limitations of this "4f-core" model and 
achieved very good agreement with experiment for the 
ground state lattice constant and the magnetic moment 
of FM bulk Gd by using a "4f-band" model and LDA, 
but he did not consider a possible AFM phase. Recent 
full potential linear muffin-tin orbitals (LMTO) calcu- 
lations show that both the LDA and the GGA 
yield an AFM phase that is lower in energy than the 
FM phase. This problem is solved by employing the 
LDA+U method to treat the electron correlations 
for the 4f-elcctrons of Gd. It was demonstrated that 
the use of LDA-I-U instead of LDA yields a correct 
FM ground state and also provides 4f electron binding 
energies in good agreement with experiment. Recently 
we confirmed quantitatively the conclusions of Ref. M] 
for bulk Gd using the LDA-I-U total energy functional 
with the full-potential linearized augmented plane wave 
(FLAPW) method 

The first total energy FLAPW calculations of the mag- 
netic ordering at the Gd surface using the 4f-core 
model reported an AFM coupling of the surface layer 
with respect to the bulk, in agreement with early in- 
terpretation of experimental data |^ . The full poten- 
tial LMTO calculations using the 4f-corc model did 
not reproduce the results of Ref. and yielded FM 
coupling between the surface and bulk magnetization, in 
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agreement with the most recent experiments It is 

thus clear that: (i) the 4f-band model with LDA fails to 
account for the correct magnetic ordering for bulk Gd, 
and the use of GGA instead of LDA does not improve 
the situation; (ii) two LDA 4f-core model calculations 
with LDA yield conflicting results ( Q and 0) for the 
magnetic ground state at the Gd surface. 

Since the LDA+U method works well to describe the 
electronic structure and magnetic ground state for bulk 
Gd, we decided to apply it in electronic structure calcu- 
lations and total energy determinations of the magnetic 
behavior of the Gd surface. 

2. Computational results. The FLAPW method 
is employed to perform scalar-relativistic self-consistent 
film calculations for Gd. The fully relativistic self- 
consistent version jl^ of this method is then used to 
perform the final LDA calculations. The LDA+U calcu- 
lations are based on the scalar-relativistic version of the 
FLAPW method @. The literature values ^ of the 
on-site repulsion U = 6.7 eV and exchange J — 0.7 eV 
were used in the calculations. 

For the Gd(OOOl) surface, we choose the isolated slab 
model based on 7-layer Gd film (with z-reflection symme- 
try) and in the first set of calculations use the bulk lat- 
tice constant (3.634 A) and c/a ratio (1.587) [Q. Here, 
32 special k-points pl| ] in the irreducible 1/3 part of 
the 2D BZ |23] were used, with Gaussian smearing for 
k-points weighting. The "mufRn-tin" radius values of 
Rmt = 3.2 a.u. and Rmt x K^ax = 9.6 (where, 
Kmax is the cut-off for LAPW basis set) were used. 

2 A. LDA results. The spin magnetic moments for a 
Gd film with its surface layer magnetically coupled par- 
allel (tt) and antiparallel (|t) to the FM bulk resulting 
from the scalar-relativistic LDA calculations show that 
the magnetically active f states are almost fully polar- 
ized (with the magnetic moment 6.88 in the bulk and 
6.82 /is at the surface) and induced spin polarization of 
« 0.5 /is/atom of conduction electrons (mainly d states). 
This magnetic coupling is a result of intraatomic inter- 
band exchange interaction between conduction band and 
localized f electrons as incorporated in the s-f exchange 
model 1^ and can be understood to be due to a positive 
inter-band d-f exchange coupling ||2^ . There is a slight 
decrease of the 4f magnetic moment at the surface layer 
due to an increase of minority spin 4f occupation. 

Starting from the results of scalar-relativistic calcula- 
tions, we then performed self-consistent relativistic LDA 
calculations for a Gd-film, assuming [0001] spin axis di- 
rection. The spin moments are slightly decreased for f 
states in comparison with scalar-relativistic calculations, 
due to an increase of minority spin occupation of the 
4f states. The small spin-orbit induced orbital magnetic 
moments (0.14 ^ib for the bulk and 0.33 /is for the sur- 
face atoms) are mainly due to the 4f minority spin contri- 
bution. A parallel coupling between spin and orbital mo- 
ments for 4f states is consistent with the 3rd Hund rule. 
The orbital moments from 5d states are about 0.02/iB 
per Gd atom and coupled anti-parallel to the spin mo- 



ments, again consistent with the 3rd Hund rule. The 
values of total magnetic moment (the sum of spin and 
orbital moments) are close to the values of spin moment 
from scalar-relativistic calculations. 

The total energy difference Ai?(||_||) between the two 
surface magnetic configurations IT and jj defined above 
is positive (36 meV/atom in scalar-relativistic calcula- 
tions) and does not change appreciably when spin-orbit 
coupling is included (40 meV/atom). The effect of the 
spin-orbit interaction is seen to be very small for the en- 
ergetics of Gd due to the fact that the 4f spin-majority 
band is fully occupied and the 4f spin-minority band is 
almost empty. Therefore, the spin-orbit coupling does 
not affect the calculated values of magnetic and total en- 
ergy properties of Gd and does not assist in resolving the 
limitations of LDA. 

2B. LDA + U results. The spin magnetic moments for 
a Gd film with surface layer magnetically coupled paral- 
lel to the FM bulk resulting from the scalar-relativistic 
LDA-I-U calculations are shown in Table |. There is a 
moderate enhancement of the magnetic moment of 4f 
states compared to LDA values (« 0.1 /is/atom) due 
to an upward shift of 1.5 eV of minority spin 4f states. 
There is practically no difference between surface and 
bulk f state magnetic moments. The total magnetic mo- 
ment at the surface layer is enhanced compared to the 
bulk mainly due to an increase of the d-state contribu- 
tion. 

The total energy difference Ai5m_||-) (71 meV/atom) 
is positive and of the same order of magnitude as the re- 
sult of the 4f-core model [Q. The results of the present 
LDA-I-U calculations for both bulk AE(^afm-fm) (63 
meV/atom p6| - the difference in energies between bulk 
AFM and FM spin configurations) and surface Ai5m_||-) 
(71 meV/atom) are in reasonable agreement with the re- 
sults of 4f-core model calculations (85 meV/atom for the 
bulk and 95 meV/atom for the surface) n^. It shows 
that parallel coupling between surface and bulk magne- 
tization is energetically preferable and there is no anti- 
parallel surface-to-bulk magnetic coupling for the Gd sur- 
face. This conclusion is consistent with experimental 
observations of the in-plane component of surface layer 
magnetization to be parallel to the bulk [||-|8| . 

The electron density of states (DOS) for the case of (en- 
ergetically preferred) || coupled surface layer are shown 
in Fig. 1. There is a 4.5 eV downward shift of the major- 
ity spin 4f states and a 1.5 eV upward shift of minority 
spin 4f states compared to the LDA calculation results. 
This latter shift makes the minority spin 4f band practi- 
cally empty and corrects the fundamental error of LDA 
4f-band model. The exchange splitting of 4f states is 
enhanced in LDA-I-U calculations by the amount of the 
Hubbard U |^6| resulting in an 11 eV splitting of ma- 
jority and minority 4f states, close to the experimentally 
derived value (12 eV) ||^. The formation of a surface state 
at the Gd surface clearly shows up as a peak of DOS in 
the vicinity of Fermi level (cf. Fig. 1) due to major- 
ity d states. From the DOS it is clear that LDA-I-U show 



2 



strongly localized character of 4f-states for both bulk and 
surface. However, the response of the 4f states to their 
environment does not allow them to be considered as a 
true core states. 

In order to check numerical convergence of our re- 
sults with respect to k-space integration, we increased 
the number of special k-points in the irreducible part of 
2D BZ from 32 to 50 in self-consistent calculations 
and found very little change in magnetic moment for both 
J,t and IT magnetic configurations (< 0.04^^). The cal- 
culated total energy difference A£'m_||-) = 72 meV per 
surface atom agrees very well with its value of 71 meV 
for a smaller number of k-points. 

3. Strain Relaxation and Magnetic Ordering at the Gd 
Surface. LEED measurements [psf show that there is 
atomic structural relaxation near the Gd surface: the in- 
terlayer distance between surface and sub-surface Gd lay- 
ers is about 2.6% smaller than its bulk value and the sub- 
surface-to-bulk layer distance is about 1% bigger than 
its bulk value. We have performed LDA-I-U calculations 
with the surface and sub-surface layers (i) with inter- 
layer distances taken from the experiment and (ii) 
with interlayer relaxations taken to be half way between 
the experimental surface values and the bulk values. As 
in the case of an ideal Gd surface we have considered 
two possible magnetic configurations (ttj it) with the 
surface layer coupled parallel and antiparallel to the FM 
bulk Gd. Here, 50 special k-points in the irreducible 1/3 
part of the 2D BZ were used. 

The surface relaxation affects very little the values of 
the magnetic moment in comparison with the ideal sur- 
face: for both cases of IT and |t coupled surface layer 
there is a slight decrease in the values of the surface and 
sub-surface layers magnetic moment due to the change 
of conduction band magnetization caused by reduced in- 
terlayer distance. 

There is, on the other hand, a surprisingly large en- 
hancement of the magnetic coupling energy A£^m_||) 
(cf.. Table ||) due to the surface relaxation: the energy 
difference increases by 90 % in the comparison with the 
unrelaxed structure. 

As was already mentioned, there is considerable exper- 
imental evidence of Tc enhancement at the Gd surface. 
Since 4f-magnetic moments are well localised and interact 
due to RKKY-type of exchange interactions, the use of 
the Heisenberg-type Hamiltonian for the dependence of 
energy on spin configuration is physically justified for Gd. 
For the sake of simplicity, we neglect the long-range be- 
havior of exchange interactions for the Gd bulk and sur- 
face, assuming that significant physics can be discussed 
in terms of nearest-neighbor (NN) interactions and ne- 
glect anisotropy in exchange interaction between a Gd 
atom and its six in-plane and six interplane NN in the 
bulk and three interplane NN at the surface. The spin 
Hamiltonian is then given by: 

= — So ^ 5,; — ^ ^ Ji,i+5S'iS'i+5 (1) 
i i 5 



where, Bq is an external field, Ji^i+s is an exchange cou- 
pling constant between the spin i and its S NN {J^ in 
the bulk and J* at the surface) and Si is a spin operator. 
We then apply "molecular field" theory |26j to Eq.(l). 
It leads to different molecular fields acting on the spin 
at the surface and in the bulk, due to the different num- 
ber of the interplane NN (six in the bulk, three at the 
surface) , and the difference between bulk and surface ex- 
change coupling constants Jt and J^. In the vicinity of 
the Curie temperature, when the value of the average 
spin moment < S{T) > is small, it is possible to intro- 
duce a Landau-Ginzburg t ype model for the temperature 
dependence of < S(T) > ]27|] using the continuum limit 
of the molecular field theory. Applying the procedure of 
Ref. to the case of the hcp(OOOl) surface, we obtain 
the result that, for 

1.5 - 2^ > or ^ > 4/3 (2) 

there is in addition to the bulk Curie temperature T^! a 
surface Curie temperature which is connected to the 
bulk Tc as: 

T! ^ [1 + (1.5 - 2:^r] Tl ■ Tl = ^'-^"^f + (3) 
Js 6kB 

The ratio A£:(||_||)/A£:(Ai^Af - FM) (cf.. Table ||) 
is then used to determine / j'' in Eq.(2) In the 

case of an ideal surface the condition Eq.(2) is not sat- 
isfied and there is no additional T^. However, when the 
surface relaxation is taken into account, the condition 
Eq.(2) is satisfied and Eq.(3) yields = 1.33 T^! in 
very good quantitative agreement with the recent exper- 
imental data § {T^ « 1.29 T^''). 

To summarize, we have presented the results of one of 
the first applications of the LDA-I-U total energy method 
to study the magnetic and electronic properties of a cor- 
related metal. We have found that the use of LDA-I-U 
instead of LDA yields FM alignment between surface and 
bulk magnetic moments, in agreement with experiment. 
An interlayer surface-to-bulk effective exchange coupling 
is calculated to be close to its bulk value for an ideal 
surface, but is enhanced by 90 % by surface relaxation. 
This enhancement is sufficiently strong to produce an 
elevated Curie temperature at the surface, as observed 
experimentally. These results also have important impli- 
cations for other rare earth (and possible also transition 
metal) surfaces and interfaces, and suggest that both 4f 
(and 3d) correlation and atomic structure must be accu- 
rately accounted for in order to quantitatively describe 
magnetism. 
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TABLE I. Spin magnetic moments {Ms in ^b) for a Gd 
film with surface layer coupled parallel to the FM bulk re- 
sulting from scalar-relativistic LDA-I-U calculations with ex- 
perimental lattice constants. For "muffin-tin" (MT), these 
values are integrals over a "muffin-tin" sphere of radius 3.2 



a.u. 


Ms layer 


s 


P 


d 


f 


total 


MT C 


0.016 


0.079 


0.46 


6.97 


7.536 


MT S-2 


0.019 


0.081 


0.49 


6.97 


7.576 


MT S-1 


0.011 


0.093 


0.50 


6.97 


7.588 


MT S 


0.041 


0.074 


0.67 


6.975 


7.773 


Interstitial: 


2.076 


Vacuum: 




0.088 





TABLE II. Total energy difference between two mag- 
netic configurations with the surface layer magnetically 
coupled antiparallel {[\) and parallel (It), Ai5(j^|_||) 
(meV/atom), for a Gd surface with (a) an ideal bulk 
atomic structure (b) a relaxed structure from experiment 
and (c) an "average" (between bulk and relaxed sur- 
face) values of surface-to-subsurface-to-bulk distances and 
its ratio to the total energy difference between AFM 
and FM bulk {AE{AFM - FM) = 63 meV/atom): 
AE(i-^_-^-^)/AE{AFM - FM) 





^-£^(iT-TT) 


A£(||_n)/A£(^FM - FM) 


a 


72 


1.14 


b 


135 


2.14 


c 


136 


2.16 
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